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Abstract 

Background: The community-associated methicillin-resistant 5. aureus (CA-MRSA) ST93 clone is becoming dominant 
in Australia and is clinically highly virulent. In addition, sepsis and skin infection models demonstrate that ST93 
CA-MRSA is the most virulent global clone of S. aureus tested to date. While the determinants of virulence have 
been studied in other clones of CA-MRSA, the basis for hypervirulence in ST93 CA-MRSA has not been defined. 

Results: Here, using a geographically and temporally dispersed collection of ST93 isolates we demonstrate that 
the ST93 population hyperexpresses key CA-MRSA exotoxins, in particular a-hemolysin, in comparison to other 
global clones. Gene deletion and complementation studies, and virulence comparisons in a murine skin infection 
model, showed unequivocally that increased expression of a-hemolysin is the key staphylococcal virulence 
determinant for this clone. Genome sequencing and comparative genomics of strains with divergent exotoxin 
profiles demonstrated that, like other 5. aureus clones, the quorum sensing agr system is the master regulator of 
toxin expression and virulence in ST93 CA-MRSA. However, we also identified a previously uncharacterized 
AraC/XylS family regulator (AryK) that potentiates toxin expression and virulence in 5. aureus. 

Conclusions: These data demonstrate that hyperexpression of a-hemolysin mediates enhanced virulence in 
ST93 CA-MRSA, and additional control of exotoxin production, in particular a-hemolysin, mediated by regulatory 
systems other than agr have the potential to fine-tune virulence in CA-MRSA. 
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Background 

Community-associated methicillin-resistant Staphylococcus 
aureus (CA-MRSA) is an emerging global problem with 
very similar clinical presentations across different clones, 
despite significant genetic diversity [1]. Many CA-MRSA 
strains carry lukSF-PV in the accessory genome, which en- 
codes the Panton- Valentine leukocidin (PVL), an exotoxin 
that causes neutrophil lysis [1]. Although there has been 
considerable controversy as to the role of this toxin in CA- 
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MRSA pathogenesis, some of this may be explained by a 
variable, species dependent susceptibility to PVL - human 
and rabbit neutrophils are lysed by PVL at very low con- 
centrations whilst mouse and monkey neutrophils are less 
susceptible, making the interpretation of animal model data 
difficult in some cases [2]. Additionally, the importance of 
PVL is also likely to be dependent on the site of infection. 
In the rabbit pneumonia model, PVL has been demon- 
strated to have a clear role in mediating severe lung necro- 
sis and inflammation [3]. In contrast, in skin infection, even 
in the rabbit model, its role remains less clear [4,5]. 

Notwithstanding PVL, the increased expression of other 
core genome virulence determinants also contributes sig- 
nificantly to the increased virulence of CA-MRSA strains 
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[6,7]. These include a-hemolysin (Hla) and a-type phenol 
soluble modulins (PSMs). Hla is a pore-forming exotoxin 
that lyses many cells including red cells, platelets, mono- 
cytes and endothelial cells [8]. Hla has been demonstrated 
to be an important mediator of virulence in skin infection 
and pneumonia [9,10]. The a-type PSMs have been re- 
cently characterized and they lyse neutrophils and red 
cells [11,12]. The a-type PSMs also mediate virulence in 
skin infection and septicemia and of these, PSMa3 is the 
most potent [11]. 

The study of unique, distantly related CA-MRSA clones 
that also demonstrate enhanced virulence, may provide 
insights into the emergence of the global CA-MRSA 
phenomenon, and also help define the genomic determi- 
nants of enhanced virulence. In Australia, the singleton 
ST93 CA-MRSA clone ST93 ("Queensland clone") has be- 
come dominant in the community [13], and we and others 
have demonstrated that our reference ST93 strain JKD6159 
was highly virulent and caused severe skin infection in a 
mouse model compared to other CA-MRSA strains includ- 
ing USA300 [14,15]. 

In this study we used exotoxin analysis, functional geno- 
mics and a murine infection model to investigate the rela- 
tive contribution of a-hemolysin, a-type phenol soluble 
modulins and Panton- Valentine leukocidin to the enhanced 
virulence of ST93 CA-MRSA. We show that increased 
virulence in the BALB/c mouse skin infection model is less 
dependent on a-type phenol soluble modulin or Panton- 
Valentine leukocidin production but is instead due to high- 
level expression of a-hemolysin in this clone, controlled 
predominantly by the agr system. 

Results and discussion 

The emergence of CA-MRSA is a major public health 
issue, and there is a clear need to understand the basis for 
both virulence and transmission of global clones of CA- 
MRSA. The genetically distinct CA-MRSA clone ST93-IV 
[2B] has rapidly become the dominant clone in Australia 
and its rise accounts for the increase in incidence of CA- 
MRSA as a whole in this country [13]. We, and others 
have previously shown that ST93 strain JKD6159 is the 
most virulent global clone of S. aureus in murine models 
[14,15]. To determine the mediators of virulence in this 
clone we initially studied exotoxin expression in a large 
collection of ST93 S. aureus from around Australia, and 
compared representative high and low expressing strains 
to an international selection of clones. 

Exotoxin expression in ST93 CA-MRSA strains 

Staphylococcus aureus expresses a wide range of exotoxins 
that may contribute to virulence. Because Hla, PVL and a- 
type PSMs have been found by others to be important viru- 
lence factors in CA-MRSA strains [9,11,16], we measured 
in vitro expression of these exotoxins by the wildtype ST93 



strains and non-ST93 comparator strains. The main 
isolates used in this study are described in Table 1, 
while the collection of ST93 isolates from around 
Australia used for comparative exotoxin expression is 
from a study by Coombs et al. [17] and summarized in 
Additional file 1. The comparison of expression of 
international clones to the ST93 reference strain 
JKD6159 and three additional ST93 strains selected for 
genome sequencing (see below) are shown in Figure 1, 
while the results for all 59 ST93 isolates compared to 
USA300 are shown in Additional file 2 (a-type PSMs) 
and Additional file 3 (Hla). The results of PVL analysis 
for the ST93 collection has been previously reported 
[17]. Because PVL is a 2-component exotoxin and both 
LukS-PV and LukF-PV are required for activity, we 
chose to measure LukF-PV expression by quantitative 
Western blot. LukF-PV was chosen over LukS-PV to 
obtain anti-LukF-PV antibody with increased specifi- 
city of binding as there was more sequence divergence 
between lukF-PV and the orthologous 2-component S. 
aureus exotoxins compared to lukS-PV. Although there 
are four a-type PSMs, PSMa3 causes the most sig- 
nificant neutrophil lysis [11] and we measured de- 
formylated and N-formylated PSMa3 expression by 
high performance liquid chromatography (HPLC). The 
proportion of deformylated and formylated forms of 
PSMs depends on the growth conditions, the activity 
of staphylococcal peptide deformylase and strain back- 
ground [7,18]. 

PVL 

As previously reported [17], PVL expression was consist- 
ent across most ST93 strains. We found that there was 
no significant difference in the LukF-PV expression in 
the PVL positive strains JKD6159, TPS3104, USA300 
and JKD6177. Although USA300 appeared to produce 
less LukF-PV than JKD6159, the difference was not sta- 
tistically significant (p = 0.0943, Figure 1A). 

PSMa3 

We found that the deformylated form of PSMa3 was al- 
most always more abundant than the N-formylated form 
(Figure IB and Additional file 2). The ST30 CA-MRSA 
strain JKD6177 did not produce any PSMa3. There was 
no significant difference in PSMa3 expression between 
JKD6159 compared to USA300, however JKD6159 pro- 
duced more PSMa3 compared to JKD6272 (p = 0.0003) 
and JKD6009 (p = 0.0003). Compared to the other ST93 
MRS A strains, JKD6159 produced more PSMa3 com- 
pared to TPS3105 (p < 0.0001), and TPS3106 (p = 0.01) 
but less than TPS3104 (p = 0.0029) (Figure IB). Expres- 
sion levels across the whole ST93 collection were vari- 
able, although many isolates produced levels at least 
equivalent to USA300 (Additional file 2). 



Table 1 Bacterial strains used in this study 


S. aureus strain 


Molecular type 


Date of isolation 


Place of isolation 


Site of isolation 


Relevant characteristics 


lukSF-PV 


reference 


Clinical isolates 
















JKD6159 


ST93-IV [2B] 


2004 


Victoria, Australia 


Blood 


Dominant Australian CA-MRSA clone 


+ 


[14] 


TPS3104 


ST93-IV [2B] 


2009 


Western Australia, Australia 


Nasal cavity 


Dominant Australian CA-MRSA clone 


+ 


This study 


TPS3105 


ST93-IV [2B] 


2005 


New South Wales, Australia 


Blood 


Australian CA-MRSA clone 




This study 


TPS3106 


ST93-V [5C2&5] 


2008 


Western Australia, Australia 


Nasal cavity 


Australian CA-MRSA clone 


- 


This study 


JI\UOZ/Z 


CT1 |\/ roDi 
j 1 1 — 1 V [ZD] 


ZUUZ 


Victoria, Australia 


DlOOQ 


Australian L-M-ivirOM cione 




L I4J 


JKD6260 


ST1-IV [2B] 


2008 


Western Australia, Australia 


Skin 


Australian CA-MRSA clone 


+ 


[14] 


JKD6177 


ST30-IV [2B] 


2003 


Melbourne, Australia 


Blood 


Australian CA-MRSA clone 


+ 


[14] 


FPR3757 USA300 


ST8-IV [2B] 


NA 


San Francisco, USA 


Wrist abscess 


Dominant North American CA-MRSA clone 


+ 


[19] 


JKD6009 


ST239-II] [3 A] 


2002 


New Zealand 


Wound 


Dominant Australian hospital-associated MRSA clone, AUS2/3 


- 


[20] 


Mutant strains 
















JrxDO 1 DylMUKjr-rV 


l\/ nRl 

o I yo-iv [Zd\ 








isogenic unmarKea iuKjr-rv wu ot jixudi jy 




This study 


JKD6159A/-7/6/ 


ST93-IV [2B] 








Isogenic unmarked hla KO of JKD6159. Deletion encompassed 
genome coordinates 1 1 21 291 -1 1 20441 . 


+ 


This study 


JKD6159Ah/a r 


ST93-IV [2B] 








Isogenic unmarked hla KO repaired in JKD6159Ah/a. Introduction 
of a novel Pstl site within hla. 


+ 


This study 


JKD6159Apsma 


ST93-IV [2B] 








Isogenic unmarked psm-a KO in JKD6159. Deletion encompassed 
genome coordinates 453364-45378. 


+ 


This study 


JKD6159Apsma r 


ST93-IV [2B] 








Isogenic unmarked psm-a KO repaired in JKD6159Apsm-a. Introduction 
of a novel Sail site within psm-a 


+ 


This study 


JKD6159A00043 


ST93-IV [2B] 








Isogenic unmarked SAA6159_00043 KO of JKD6159. Deletion encompassed 
genome coordinates 53156 - 54561 


+ 


This study 


JKD6159_AraCr 


ST93-IV [2B] 








Isogenic AraC/XylS regulator repaired in JKD6159 


+ 


This study 


TPS3105r 


ST93-IV [2B] 








Isogenic agrA repair of TPS3105 




This study 



KO: knockout, NA: not available. 
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Figure 1 /n vitro exotoxin expression of wildtype S. aureus 
isolates. JKD6159 (ST93-IV [2B]) compared with non-ST93 CA-MRSA 
strains FPR 3757 USA300 (ST8-IV [2B]), JKD6177 (ST30-IV [2B]) ; and 
JKD6272 (ST1-IV [2B]); Hospital-associated MRSA strain JKD6009 
(ST239-III [3A]) ; wildtype ST93 strains TPS3104 (ST93-IV [2B]), TPS3105 
(ST93-IV [2B]) ; and TPS3106 (ST93-V [5C2&5]). (A) LukF-PV expression 
measured by quantitative Western blot. RN4220 was included as a 
negative control because it does not contain lukF-PV. All PVL negative 
strains did not express LukF-PV. There was no significant difference in 
the amount of LukF-PV expressed by the S. aureus strains containing 
lukSF-PV. Data shown are mean intensity of bands in arbitrary units 
and SEM. (B) PSMa3 expression measured by HPLC. JKD6177 did not 
produce PSMa3. JKD6272 (p = 0.0003), JKD6009 (p = 0.0003), TPS31 05 
(p < 0.0001) and TPS31 06 (p = 0.01 00) produced less deformylated and 
N-formylated PSMa3 compared to JKD61 59. There was no difference 
between PSMa3 production byJKD6159and USA300. TPS3104 
expressed more PSMa3 than JKD6159 (p = 0.0029). Data shown are 
mean concentration (ug/ml), presented as vertical stacked bars and 
SEM. Deformylated PSMa3 is shown in grey bars. N-formylated PSMa3 
is shown in white bars. (C) Hla expression measured by quantitative 
Western blot. RN4220 was included as a negative control because it 
does not express Hla. JKD6159 expressed more Hla compared to all 
non-ST93 wildtype strains (p < 0.0001 for all strains except JKD61 77 
p = 0.0107). TPS31 05 and TPS3106 produced significantly less Hla 
(p < 0.0001). There was no difference in Hla production between 
JKD6159 and TPS3104. Data shown are mean intensity of bands in 
arbitrary units and SEM. Note, ***p < 0.001, **p < 0.01, *p < 0.05. 



Hla 

Hla expression appeared high for the majority of ST93 
isolates, with the exception of four strains where expres- 
sion was low (Additional file 3). JKD6159 produced 
greater levels of Hla than all the wildtype strains, inclu- 
ding USA300 (p < 0.0001 for all strains except JKD6177, 
p = 0.0107, Figure 1C). There was no difference in Hla 
expression between JKD6159 and TPS3104. 

Here we have demonstrated that the majority of ST93 
strains consistently produce higher levels of Hla com- 
pared to other clones, including USA300, while produc- 
tion of PVL and a-type PSM is similar, suggesting that 
enhanced expression of Hla may be responsible for in- 
creased virulence of ST93 CA-MRSA. 



Comparative virulence of ST93 isolates with differential 
exotoxin expression 

To further examine the role of these selected exotoxins in 
our mouse skin infection assay, we compared the virulence 
of four ST93 isolates selected based on their exotoxin 
expression profiles (high exotoxin expression, JKD6159 
and TPS3104; low exotoxin expression, TPS3105 and 
TPS3106). TPS3104 was as virulent as JKD6159 in the 
mouse model in all outcome measures (Figure 2). In 
contrast, the strains with reduced exotoxin expression 
TPS3105 and TPS3106 were significantly less virulent 
compared to JKD6159, with less weight loss at day 5 of in- 
fection (p < 0.0001), smaller lesion size (p < 0.0001) and 
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Figure 2 Virulence characteristics of wildtype ST93 CA-MRSA isolates. S. aureus JKD6159 compared with three other wildtype ST93 
CA-MRSA isolates, TPS3104, TPS3105 and TPS3106 in a BALB/c mouse skin infection assay. At least 10 mice were used for each bacterial strain. 
(A) Weight loss induced by intradermal infection with S. aureus strains is demonstrated as percentage loss of weight over 5 days. The difference 
in percentage weight loss between JKD6159 and TPS3105 and TPS3106 was significant (p < 0.0001). There was no difference in weight loss 
between JKD6159 and TPS3104. Data shown are mean weight loss and SEM. (B) Skin lesion area (mm 2 ) at 5 days after infection was significantly 
greater with JKD6159 infected mice compared to TPS3105 and TPS3106 (p < 0.0001). There was no difference in lesion area between JKD6159 
and TPS3104. Data shown are mean area and SEM. (C) Recovery of S. aureus (log CFU) from infected tissues at 5 days after infection from 
JKD6159 infected mice was greater than with TPS3105 (p = 0.01 77) and TPS3106 infected mice (p = 0.0328). There was no difference between 
JD6159 and TPS3104 infected mice. Data shown are mean CFU and SEM. Note, ***p < 0.001, *p < 0.05. 



less CFU recovery from lesions (TPS3105, p = 0.0177; 
TPS3106, p = 0.0328) in the model (Figure 2). 

Impact of exotoxin expression on virulence of ST93 
CA-MRSA in the murine skin infection model 

To further characterize the contribution of each of the 
exotoxins to disease in the murine model, genetic deletion 
and complementation experiments were performed for 
each of the selected toxins. 

Hla 

Given the increased in vitro expression of Hla by JKD6159 
and TPS3104 and the apparent correlation of this increased 
expression with increased virulence in the mouse skin in- 
fection model, we generated ]KD6159Ahla and assessed 
this mutant in the mouse skin infection assay (Figure 3). 
There was a marked attenuation in virulence in all out- 
come measures with significantly decreased weight loss 
(p < 0.0001), lesion size (p < 0.0001) and CFU recovery 
(p = 0.0177). To confirm that an unintentional mutation 
introduced during the procedure to knock-out hla was 
not responsible for the reduced virulence in this strain, 
complete genome sequencing of the strain using Ion 
Torrent sequencing was performed. Mapping of se- 
quence reads from JKD6159AW<z against JKD6159 (40x 
genome coverage) demonstrated no additional differences 
between JKD6159 and JKD6159AW*. We also restored Hla 
function in this strain by using allelic replacement to repair 
the deletion in hla and create strain ]KD6159Ahla r. As ex- 
pected, Hla expression was absent in ]KD6159Ahla and ex- 
pression was restored in ]KD6159Ahla r when tested by 
Western Blot (Additional file 4A). JKD6159A/zk r also 
reverted to high virulence in the mouse skin infection assay 
(Figure 3). The apparent slight reduction in virulence of 
this hla repaired strain compared to wild type JKD6159 is 



explained by incomplete penetration of the restored hla 
allele in ]KD6159Ahla r, resulting in mixed bacterial popu- 
lations and reversion to ]KD6159Ahla for some of the 
mice (Additional file 4B and C). 

a-type PSMs 

In order to determine the contribution of a-type PSMs to 
virulence of JKD6159, we generated JKD6159 Apsma 
(deletion of the whole a-type PSM locus) and assessed 
this mutant in the mouse skin infection assay (Figure 3). 
There was no significant difference in virulence in all 
outcome measures; weight loss (p = 0.06), lesion size 
(p = 0.8174) and CFU recovery (p = 0.1925). The absence 
of PSMa3 production was verified by HPLC and the in- 
tegrity of the JKD6159 Apsma was confirmed by com- 
plementation (Additional file 5). It is noteworthy that 
the level of expression of PSMa3 by JKD6159 was simi- 
lar to USA300 (Figure 1), a strain that produces high 
levels of PSMs and where a contribution to virulence 
has been demonstrated [7,11]. Despite this, the deletion 
mutant (JKD6159 Apsma) demonstrated no attenuation 
of virulence compared to JKD6159 (Figure 3). The sig- 
nificantly divergent genetic background of ST93 com- 
pared with USA300 may account for this difference 
in the importance of a-type PSMs to the virulence of 
JKD6159 [6]. 

PVL 

We constructed an isogenic PVL negative mutant in 
JKD6159 by deleting lukSF-PV. Western Blot analysis con- 
firmed the absence of LukF-PV in the mutant (Additional 
file 6). Assessment of the JKD6159 AlukSF-PV mutant 
in the mouse skin infection model showed no decrease 
in virulence (Figure 3). Therefore PVL was not contrib- 
uting to the increased virulence in JKD6159 in this 
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A B C 




Figure 3 Virulence characteristics of S. aureus JKD6159 and isogenic exotoxin mutants derived from JKD6159. JKD6159 compared 
to isogenic PVL knockout (JKD61 59 MukSF-PV), isogenic Hla knockout (JKD61 59 Ahla), isogenic Hla complemented strain (JKD61 59Ahla r) and 
isogenic PSM-a knockout (JKD6159Apsma) in a BALB/c mouse skin infection assay. (A) Weight loss induced by intradermal infection with 
S. aureus strains is demonstrated as percentage loss of weight over 5 days. There was no significant difference between JKD6159, 
JKD6159A/L//CSF-P1/ and JKD6159Apsma infected mice. There was significantly less weight loss in mice infected with JKD61 59Ahla compared 
to JKD6159 (p< 0.0001). There was also less weight loss in mice infected with JKD6159Ah/a compared to JKD6159A/i/a r (p = 0.0063). Mice 
infected with JKD6159A/i/c/ r had less weight loss compared to JKD6159 (p = 0.0004). Data shown are mean weight loss and SEM. (B) 
There was no difference in skin lesion area (mm 2 ) at 5 days after infection in mice infected with JKD6159 and JKD6159A/u/c5/ r -/V and 
JKD6159Apsma. Mice infected with JKD6159A/i/a had significantly smaller lesions (p < 0.0001). In some mice, there was no cutaneous lesion 
seen. There were significantly smaller lesions in mice infected with JKD6159A/i/c? compared to JKD6159A/i/c? r (p < 0.0001). Mice infected with 
JKD6159A/i/6/ r had smaller lesions compared to JKD6159 (p = 0.024). Data shown are mean area and SEM. (C) Recovery of S. aureus (log CFU) 
from infected tissues at 5 days after infection from JKD6159 infected mice was no different to that from JKD6159A/i//c5F-Pl/, JKD6159Ap5/na and 
JKD6159A/7/C/ r. There was significantly less S. aureus recovered from JKD6159A/i/a infected mice (p = 0.0177). There was also significantly less S. aureus 
recovered from JKD6159A/i/a infected mice compared to JKD6159A/i/a r (p = 0.0018). Data shown are mean CFU and SEM. Note, ***p < 0.001, *p < 0.05, 
compared to JKD6159. 



murine model Murine neutrophils, unlike rabbit and 
human neutrophils are relatively resistant to the effects 
of PVL so it is difficult to draw firm conclusions as to 
the human importance of this result [2]. However, the 
aim of this study was to uncover the mechanisms for 
the observed increased virulence of ST93 previously 
demonstrated using this mouse model [14]. Our results 
reinforce the results of others who have examined dif- 
ferent S. aureus clones which indicate that Hla, rather 
than PVL is the main mediator of virulence in CA-MRSA 
in a mouse skin infection model [9,10,21,22]. It should be 
noted that other authors have concluded that the rabbit 
skin infection model gave very similar results to the mouse 
model for infection at the same site [4]. Nonetheless, testing 
of our PVL deletion mutant in a rabbit model may be 
warranted in future. 

Genome sequencing of three additional ST93 isolates 

We have previously fully sequenced and annotated the 
genome of ST93 strain JKD6159 [14,23]. The differential 
virulence and exotoxin expression of some ST93 isolates 
compared to JKD6159 was then exploited by using whole 
genome sequencing and comparative genomics to deter- 
mine the genetic basis for exotoxin expression in this 
clone. We selected the high expression strain TPS3104 
and the low virulence and expression strains TPS3105 and 
TPS3106 to compare to JKD6159. De novo assembly of 
each of these strains resulted in -700 contigs per isolate, 



with a genome length of 2.8 Mbp. The de novo assembly 
metrics are summarized in Additional file 7. The contigs 
were aligned to JKD6159 using BLASTN, with some im- 
portant differences demonstrated between the strains 
(Figure 4A). TPS3104 contained SCCmecW and <|>SA2 
with lukSF-PV; TPS3105 contained SCCmecW but lacked 
(|)SA2 and lukSF-PV; TPS3106 contained SCCmecV, and 
())SA2 without lukSF-PV. In addition, read-mapping 
against the complete JKD6159 genome (chromosome and 
plasmid: 2,832,164 bp) was then employed to define a 
2,720,685 bp core genome among these four isolates (96% 
of the JKD6159 genome), revealing 253 polymorphic nu- 
cleotide positions, some of which were common to two 
strains. A phylogeny was inferred that confirmed the close 
relationship among all isolates, with TPS3106 more dis- 
tantly related to the others (Figure 4B). The unmapped 
reads from each isolate were also subjected to de novo as- 
sembly to identify DNA not present in JKD6159. TPS3104 
and TPS3105 contained no new sequences, while TPS3106 
contained 34 kb of additional DNA, predominantly span- 
ning the SCCmecV region. 

Comparative genomics of ST93 and the importance of agr 
in the virulence of ST93 CA-MRSA 

We next explored the contribution of specific mutations to 
the differential virulence of the ST93 strains. Using our 
read mapping approach described above, we compared the 
genome sequences of TPS3104, TPS3105 and TPS3106 
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l 



OSA2 

Figure 4 Whole genome sequence analysis and comparison of JKD6159 with other ST93 CA-MRSA isolates. (A) Circular diagram of the 
JKD6159, TPS3104,TPS3105 and TPS3106 chromosomes (from inner to outer circles). TPS3 104, TPS3 105 and TPS3106 contigs were mapped by 
BLASTN to JKD6159. TPS3104 contained SCCmecIV and 4>SA2 with lukSF-PV; TPS3105 contained SCCmecIV but lacked <|>SA2 and lukSF-PV; TPS3106 
contained SCCmecV, and cJ)SA2 without lukSF-PV. (B) ST93 S. aureus phylogeny inferred by split decomposition analysis from pairwise comparisons 
of the 253 variable nucleotide positions identified from the ST93 core chromosome of 2,720,685 bp. Figures indicate the number of nucleotide 
substitutions per branch. All nodes have 100% bootstrap support. 



with each other and with JKD6159. There were a number 
of single nucleotide polymorphisms (SNPs) and insertions 
and deletions (indels) differentiating the strains from 
JKD6159 (Additional files 8, 9, 10). We searched for muta- 
tions in regulatory genes that could potentially explain the 
different virulence phenotypes of the strains. Notably, both 
avirulent ST93 strains, TPS3105 and TPS3106 contained 
mutations within the agr locus. We have since completed 
whole genome sequencing of TPS3151 and TPS3161 and 
found they contain predicted amino acid substitutions in 
AgrC that might disrupt agr function (Stinear et al, sub- 
mitted). These isolates demonstrated low expression of 
Hla (Additional file 3). 

Additionally, TPS3106 also contained a mutation in 
a gene encoding a previously uncharacterized AraC/ 
XylS family regulatory protein. This was also of par- 
ticular interest as members of this class have been 
shown to contribute to the regulation of exotoxin ex- 
pression [24,25]. 

TPS3105 contained a frame-shift mutation within agr A 
(Sa_JKD6159 nucleotide 2096502) and a further substi- 
tution (G to A) within agrA at nucleotide 2096569), 
while TPS3106 contained an -356 bp deletion spanning 
the agr effector molecule, RNAIII (deletion spanning nu- 
cleotides 2093372 to 2093728). These mutations sug- 
gested these isolates were agr deficient. To demonstrate 
that these changes explained the reduced exotoxin ex- 
pression and loss of virulence in TPS3105 we repaired 



agrA using allelic exchange to create TPS3105r (the par- 
ental strain TPS3105 now containing an intact agrA). 
Quantitative real time RT-PCR for RNAIII demonstrated 
that TPS3105r produced 325-fold more RNAIII than 
TPS3105. Virulence was also restored and TPS3105r 
caused greater weight loss, skin lesion area and CFU re- 
covery from lesions compared to the parental strain 
TPS3105 (p < 0.0001, Figure 5). There was no significant 
difference between JKD6159 and TPS3105r in all outcome 
measures in the mouse skin infection model (Figure 5). 
These experiments show that intact agr is essential for 
the virulence of ST93 CA-MRSA. The agrA repaired 
mutant of TPS3105, TPS3105r expressed significantly 
greater amounts of PSMa3 (p < 0.0001) and Hla (p = 
0.0019), consistent with agr control of these virulence 
determinants (Figure 6). Thus, despite the genetic diver- 
gence of ST93 from other S. aureus [14], the molecular 
foundation of virulence for this CA-MRSA clone is 
similar in this respect to USA300 [9,26,27] and other S. 
aureus strains [28,29], where the importance of agr has 
been very well established. 

The AraC/XylS regulator (AryK) enhanced Hla expression 
and virulence in ST93 CA-MRSA 

The SNP at position 92551 in SAA6159_00084 introduced 
a premature stop codon and created a pseudogene within 
SAA6159_00084 in JDK6159, however the gene was intact 
in TPS3106. The intact version of this gene, which was also 
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Figure 5 The importance of agr and aryK in the virulence of ST93 CA-MRSA. Isogenic repaired agr mutant TPS31 05r compared to TPS31 05 
and JKD6159, and JKD6159 compared with isogenic repaired AraC/XylS family regulator mutant (JKD61 59_AraCr) in a BALB/c mouse skin 
infection assay. At least 10 mice were used for each bacterial strain. (A) Weight loss induced by intradermal infection with S. aureus strains is 
demonstrated as percentage loss of weight over 5 days. There was no difference between JKD6159 and TPS31 05r in all outcome measures. 
TPS3105r infected mice had significantly increased weight loss compared to TPS3105 (p < 0.0001). There was a small increase in weight loss in mice 
infected with JKD6159_AraCr compared to JKD6159 (p = 0.0311). Data shown are mean weight loss and SEM. (B) Skin lesion area (mm 2 ) at 5 days after 
infection in TPS3105r infected mice was significantly increased compared to TPS3105 (p < 0.0001). Mice infected with JKD6159_AraCr had increased 
lesion area compared with JKD6159 (p < 0.0001). Data shown are mean area and SEM. (C) Recovery of S. aureus (log CFU) from infected tissues at 
5 days after infection from TPS31 05r was significantly greater than from TPS 3105 infected mice (p < 0.0001). There was no difference in S. aureus 
recovered from mice infected with JKD6159 and JKD6159_AraCr. Data shown are mean CFU and SEM. Note, *** p < 0.001, * p < 0.05. 



intact in 19 other publically available S. aureus genome 
sequences we examined, encodes a previously uncharac- 
terized AraC/XylS family regulatory protein. While the 
virulence attenuation in TPS3106 was likely a direct re- 
sult of the agr deficiency, we also wanted to determine 
if the novel regulator mutation in SAA6159_00084 im- 
pacted the virulence in ST93 S. aureus. 
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Figure 6 In vitro PSMa3 and Hla expression of mutant S. 
aureus isolates. JKD6159 compared with JKD6159_AraCr. TPS3105 
compared with TPS3105r. (A) PSMa3 expression measured by HPLC. 
JKD6159_AraCr expressed more PSMa3 than JKD6159 (p = 0.0325). 
TPS3105r expressed more PSMa3 thanTPS3105 (p< 0.0001). Data 
shown are mean concentration (ug/ml), presented as vertical 
stacked bars and SEM. Deformylated PSMa3 is shown in grey bars. 
N-formylated PSMa3 is shown in white bars. (B) Hla expression 
measured by quantitative Western blot. There was a small but 
statistically significant increase in Hla production by JKD6159_AraCr 
(p = 0.0473). TPS3105r expressed more Hla than TPS3105 (p = 0.001 9) 
Data shown are mean intensity of bands in arbitrary units and SEM. 
Note, *p < 0.05, compared to JKD6159. Note also, ### p < 0.001 and 
## p < 0.01, compared to TPS3105. 



To test the hypothesis that SAA6159_00084 encoded a 
regulator of virulence, we repaired the premature stop 
codon in SAA6159_00084 in JKD6159 using allelic ex- 
change to generate strain JKD6159_AraCr. To confirm we 
had not introduced additional DNA changes during al- 
lelic exchange we sequenced the whole genome of 
JKD6159_AraCr and found no additional mutations (35 x 
coverage). JKD6159_AraCr encoding an intact copy of 
SAA6159_00084 demonstrated a modest, but significant 
increase in virulence as indicated by lesion size (p < 
0.0001) and weight loss in the mouse skin infection 
assay (p = 0.0311, Figure 5), suggesting that this pro- 
tein is a positive regulator of virulence in CA-MRSA 
strains. JKD6159_AraCr expressed more PSMa3 (p = 
0.0325) and Hla (p = 0.0473) than its parental strain 
JKD6159 that was consistent with an increase mouse 
skin lesion size (Figure 6). We propose the name aryK 
for SAA6159_00084 (AraC family-like gene). 

RNAseq demonstrates global regulatory impact of AryK 

To investigate the regulatory impact of AryK, RNAseq 
was performed using RNA extracted from stationary 
phase cultures (Figure 7). This growth phase was se- 
lected as we reasoned that AryK might be interacting 
with agr and thus any impacts on Hla expression would 
be greatest at this time. A small number of virulence- 
associated loci were down regulated in the aryK mutant 
(JKD6159), including beta- type phenol soluble modulins 
(SAA6159_01024 and SAA6159_01025), and the viru- 
lence regulator saeS. However, the most dramatic and 
significant transcriptional changes were found in genes 
involved in central metabolic functions. Using the Kyoto 
Encyclopedia of Genes and Genomes (KEGG) pathway 
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Figure 7 Heatmap of RNA Sequencing comparing JKD6159 {aryK inactive) to JKD6159_AraCr {aryK intact). RNA seq was performed in 
duplicate from stationary phase cultures. This heatmap, clustered on expression profiles, was created based on log 2 transformed counts to 
identify consistent changes in expression profiles between strains. To be included in the heat map, genes were required to have at least 1000 
counts (reads), totaled over all samples, where the standard deviation of log 2 expression differences had to exceed two. The heatmap highlights 
significant 6/ry/(-dependent changes, in particular genes involved in the regulation of central metabolic functions. 



analysis (www.genome.jp) [30] the major pathways where 
differential transcription occurred were down regulation 
of genes linked to purine metabolism (purK, purS, purQ, 
purL, purF), valine, leucine and isoleucine biosynthesis 
(leuA, leuB, leuC), and oligopeptide transport (opp3C, 
opp3D, opp3F, opp3A). In contrast up regulation of genes 
encoding cation transport systems (rnnhB_l, mnhC_l, 
mnhD_l, mnhF_l, mnhG_l) was found. 

Here, we have clearly demonstrated that agr is the major 
"on-off switch for virulence in ST93 CA-MRSA, but we 



also found that other genetic changes are impacting 
virulence gene regulation in a clone-specific manner. 
We speculate that the inactivation of aryK may have 
been an evolutionary response by ST93 CA-MRSA to 
modulate or fine-tune the amount of Hla and other fac- 
tors required for host persistence. There are six AraC/ 
XylS family regulators in S. aureus (SA0097, SA0215, 
SA0622, SA1337, SA2092, SA2169; S. aureus strain N315 
locus tags). Two of these, Rbf (SA0622) and Rsp (SA2169) 
have been studied and demonstrated in other S. aureus 
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strains to regulate biofilm formation and modulate expres- 
sion of surface-associated proteins [24,25,31]- In contrast, 
we found that aryK increases Hla expression and virulence, 
acting as a positive regulator of virulence by directly or in- 
directly upregulating exotoxin expression, without an ap- 
parent effect on agr expression in stationary phase. 

Conclusions 

In this study, we have obtained insights into the genetic 
basis for the increased virulence of ST93 by using a 
combination of comparative and functional genomics. 
We have demonstrated the key role of Hla and agr and 
shown how an additional novel regulatory gene, aryK 
by a loss-of-function point mutation, is modulating 
virulence in this clone. Quantification of exotoxin ex- 
pression in a larger collection of ST93 strains demon- 
strated that the findings in strain JKD6159 are relevant 
to the majority of the ST93 population isolated from 
around Australia as exotoxin expression in JKD6159 is 
representative of most of the ST93 population. Our 
study highlights the power of comparative genomics to 
uncover new regulators of virulence but it also shows 
the complex nature of these changes even in closely re- 
lated bacterial populations. Careful strain selection, de- 
tailed comparative genomics analyses, and functional 
genomic studies by creating multiple genetic changes in 
one strain will be required to gain a full insight into the 
genetic basis for the emergence and hypervirulence of 
ST93 CA-MRSA. 

Methods 

Ethics statement 

This study was performed in accordance with the Australian 
Prevention of Cruelty to Animals Act 1986 and the 
Australian code of practice for the care and use of animals 
for scientific purposes. The protocol was approved by the 
Animal Ethics Committee of the University of Melbourne 
(Permit Number: 0911248.2). 

Bacterial strains and culture 

Bacterial strains used in this study are summarized in 
Table 1 (international clone collection) and Additional file 
1 (ST93 strain collection), and include the ST93 reference 
strains JKD6159, USA300 strain FPR3757 [19], ST30 strain 
JKD6177, and the HA-MRSA ST239 clone JDK6009 [20], 
as well as 58 additional ST93 collected from around 
Australia and previously reported [17]. For all experi- 
ments except exotoxin expression bacteria were grown 
in brain heart infusion broth (BHI, Oxoid). For the 
mouse skin infection assay, S. aureus were harvested at 
the stationary phase of growth after 18 hours incuba- 
tion (OD 600 ~ 2.0), washed, diluted and resuspended in 
PBS. The bacterial inoculum (CFU) and viable counts 



were determined by plating onto BHI agar and colony 
enumeration. 

For LukF-PV expression experiments, bacteria were 
grown in CCY media (3% yeast extract (Oxoid), 2% 
Bacto Casamino Acids (Difco), 2.3% sodium pyruvate 
(Sigma- Aldrich), 0.63% Na 2 HP0 4 , 0.041% KH 2 P0 4 , pH 6.7). 
For a-hemolysin (Hla) and PSMa3 expression experi- 
ments, bacteria were grown in tryptone soy broth (TSB, 
Oxoid). Overnight cultures were diluted 1:100 into fresh 
media and then incubated at 37°C with shaking (180 rpm) 
until stationary phase was achieved. For LukF-PV detection, 
isolates were cultured for 8 hours (OD 600 ~ 1.8); for 
Hla detection, isolates were cultured for approximately 
3 hours (OD 600 ~ 1.8); and for PSMa3 detection, isolates 
were cultured for 24 hours (OD 600 ~ 2.0). Culture superna- 
tants were harvested by centrifugation and filter sterilized. 
These were performed in at least triplicate for each S. 
aureus strain. 

Detection of LukF-PV and Hla by western blotting 

Trichloroacetic acid was added to culture supernatants 
and incubated at 4°C overnight. Precipitates were then 
harvested by centrifugation, washed with acetone, air- 
dried and solubilized in a SDS and 2-mercaptoethanol 
containing sample buffer. The proteins were separated on 
12% SDS-PAGE. A peptide sequence specific to LukF-PV, 
HWIGNNYKDENRATHT was synthesized and HRP 
conjugated polyclonal chicken IgY was raised against 
this peptide (Genscript). This antibody was used to de- 
tect LukF-PV with enhanced chemiluminescence. Im- 
ages generated from the western blots were quantitated 
using GS800 Calibrated Densitometer (BioRad) and 
Image J [32]. 

Hla was detected using a polyclonal rabbit anti-Hla 
(Sigma-Aldrich), in buffer containing 20 mM DEPC to 
inhibit non-specific protein A binding and HRP conju- 
gated sheep anti-rabbit secondary antibody (Millipore) 
with enhanced chemiluminescence detection [33]. For 
comparison of JKD6159 versus international clone col- 
lection (Table 1), images generated from the Hla western 
blots were quantitated using GS800 Calibrated Densito- 
meter (BioRad) and Image J [32] . 

Subsequently, for comparison of JKD6159 and other 
ST93 strains (Table 1), detection of chemiluminescence 
was performed using the MF-ChemiBIS 3.2 platform 
(DNR Bioimaging systems). Quantitation was performed 
using Image J [32]. 

Detection of PSMa3 expression 

HPLC chromatography was performed on an Agilent 
Technology 1200 Series system with an analytical Agilent 
Eclipse XDB-C18 (4.6 mm x 150 mm) column. A water/ 
acetonitrile gradient (0.1% trifluoroacetic acid) from 0 - 
100% acetonitrile over 28 min at a flow rate of 1 mL/min 
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was used. The total run time was 32 min, and peaks were 
quantified at a wavelength of 214 nm. The deformylated 
and formylated form of PSMa3 MEFVAKLFKFFKDLL 
GKFLGNN was identified in the S. aureus TSB culture su- 
pernatants by comparison of their retention times to a 
commercially synthesized PSMa3 standard (GenScript) 
and by spiking the samples with the synthesized standards. 
The identity of the deformylated peptide present in the 
samples was confirmed by analysing collected fractions by 
ESI-MS. There was only one peptide present in this frac- 
tion; the deformylated form of PSMa3. In contrast, other 
peptides were observed in the fractions of USA300, 
JKD6272, TPS3104, TPS3105r, and JKD6159_AraCr con- 
taining the N-formylated form of PSMa3. In these cases, 
the percentage of N-formylated PSMa3 peptide was deter- 
mined using the total ion count of the major peaks in the 
ESI-MS and the peak area of the HPLC chromatogram 
was adjusted accordingly. The concentrations of the defor- 
mylated and formylated forms of PSMa3 were determined 
by comparison of their peak areas to those of the synthe- 
sized standards. The standard curves were constructed in 
the concentration range of 6.2 - 100 \ig/ mL and were lin- 
ear over this range. 

DNA methods, molecular techniques and construction 
of mutants 

DNA was extracted using the GenElute kit according 
to the manufacturers instructions (Sigma- Aldrich). A 
lukSF-PV knockout, hla knockout and a repaired agrA 
of TPS3105 were generated according to the published 
method [34]. For the knockouts, flanking sequences 
were amplified and ligated prior to cloning with pKORl. 
For allelic replacement to generate TPS3105r, a PCR 
product of agrA from JKD6159 was cloned with pKORl. 
For allelic replacement JKD6159_AraCr, a PCR product 
of this AraC regulator from TPS3106 was cloned with 
pKORl. The deletion of the whole psma locus in JKD6159, 
chromosomal restoration of psma in ]KD6159Apsma and 
the restoration of Hla expression in ]KD6159Ahla were 
conducted using the pIMAY protocol described by 
Monk et al. [35]. Knockout and restoration amplimers 
were cloned into pIMAY by SLIC [36]. The primers 
used are listed in Additional file 11. The knockout and 
restoration clones were confirmed by PCR and Sanger 
sequencing of the mutated locus. Mutations were fur- 
ther validated with functional assays of activity, which 
included sheep blood hemolysis, western blot, and/or 
HPLC. 

Mouse skin infection assay 

Mice were infected with S. aureus as previously de- 
scribed [14]. Briefly, six- week-old female BALB/c mice 
were infected by intradermal injection with 10 8 CFU of 
S. aureus. Mice were assessed and weighed daily for five 



days. Mice were culled on the 5th day and lesion size 
measured and CFU recovered from infected tissues by 
homogenization and colony enumeration on BHI. For 
each S. aureus strain, at least 10 mice were assessed. 

Genome sequencing 

Genome sequences for three ST93 strains (TPS3104, 
TPS3105, TPS3106) were obtained from an Illumina 
GAIIx analyzer using 100 bp paired-end chemistry with a 
mean fold coverage of 331 x. Genome sequencing of the 
two laboratory-induced mutants JKD6159A/zk (TPS3265) 
and JKD6159_AraCr (TPS3268) was performed using Ion 
Torrent sequencing technology. 

Comparative genomics 

A read mapping approach was used to compare the se- 
quences from all isolates used in this study, as previously 
described [14,37]. Briefly, the reads from all genomes 
were aligned to the JKD6159 reference using SHRiMP 
2.0 [38]. SNPs were identified using Nesoni v0.60 [www. 
bioinformatics.net.au]. Using the whole genome sequence 
of JKD6159 as a reference, a global SNP analysis was 
performed, and allelic variability at any nucleotide position 
was tallied to generate a global SNP analysis for every gen- 
ome compared to JKD6159. 

Quantitative RT-PCR for RNAIII expression 

To investigate activity of the agr locus (RNAIII) qRT-PCR 
was performed for RNAIII as previously described [37]. 
Briefly, RNA was prepared as previously described with two 
on-column DNase I digestion steps and cDNA synthesis 
using Superscript II reverse transcriptase (Invitrogen). Rela- 
tive expression was determined as previously described and 
was normalised against gyrB. Results were obtained from 3 
biological replicates each performed in triplicate. 

RNA sequencing 

Staphylococcus aureus strains JKD6159 and JKD6159_Ar- 
aCr were grown to early stationary culture as described 
above. For RNA protection, 0.5 volumes of RNAlater® 
RNA stabilization reagent (Qiagen) was added immediately 
to the liquid culture and allowed to incubate with the bac- 
terial suspension for 15 minutes at room temperature. Cells 
were pelleted at 5,000 x g for 5 minutes followed by RNA 
extraction using RNeasy mini kit (Qiagen) and two rounds 
of DNase I digestion (Qiagen) according to the manu- 
facturer s instruction. RNA concentration was quantified 
using Qubif 2.0 Fluorometer and RNA quality assessed 
using Agilent 2100 Bioanalyzer. Ten \ig of total RNA 
from the stationary growth phase with RNA intergrity 
number (RIN) greater than 7 was used in RNA-seq. 
Ribosomal depletion, cDNA library preparation and 
pair ended sequencing using HiSeq2000 sequencing 
platform was performed by Beijing Genome Institute 
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(Hong Kong, China). RNAseq was performed on two bio- 
logical samples for each strain. 

RNAseq reads were mapped onto the JKD6159 refe- 
rence genome [23], using SHRiMP 2.2.2 [39]. Alignment 
to CDS features from each biological replicate of each 
strain provided counts that were a measure of mRNA 
levels. Counts were normalized using the trimmed-mean 
normalization function in edgeR, part of the BioConduc- 
tor package [40]. A heat map was created based on log 2 
transformed counts to identify consistent changes in ex- 
pression profiles between strains. To be included in the 
heat map, genes were required to have at least 1000 
counts, totaled over all samples, where and the standard 
deviation of the log 2 expression levels had to exceed two. 

Statistical analysis 

Percentage mouse weight change at day 5, viable counts 
of S. aureus in mouse tissues and skin lesion area of 
each isolate, Hla, LukF-PV and PSMa3 expression versus 
JKD6159 were analyzed using an unpaired t test. A simi- 
lar analysis was used to analyze virulence outcome mea- 
sures and exotoxin expression between TPS3105 and 
TPS3105r. (There was no difference in results when 
Bonferonni analysis was performed). All analyses were 
performed using Prism 5 for Macintosh v5.0b (GraphPad 
Software Inc.). 

Availability of supporting data 

The data sets supporting the results of this article are in 
the NCBI Sequence Read Archive under study accession 
SRP004474.2 and the NCBI BioProject Archive under 
study accession PRJNA217697. 

Additional files 



this strain resulted in restoration of formylated PSMa3 expression. In all 
strains 6-toxin expression was maintained. 

Additional file 6: LukF-PV Western Blot of JKD61 59 and 

JKD61 S9MukSF-PV. Western Blot demonstrating that JKD61 59NukSF-PV 

does not express LukF-PV. 

Additional file 7: Table of de novo assembly characteristics for S. 
aureus strains TPS3104, TPS3105 and TPS3106. 

Additional file 8: Table of single nucleotide differences between 
JKD6159 and TPS3104. 

Additional file 9: Table of single nucleotide differences between 
JKD6159 and TPS3105. 

Additional file 10: Table of single nucleotide differences between 
JKD6159 and TPS3106. 

Additional file 11: Table of primers used in this study. 
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